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Abstract 

We probe charge photogeneration and subsequent recombination dynamics in neat regioregular 
poly(3-hexylthiophene) films over six decades in time by means of time-resolved photoluminescence 
spectroscopy. Exciton dissociation at 10K occurs extrinsically at interfaces between molecularly 
ordered and disordered domains. Polaron pairs thus produced recombine by tunnelling with dis- 
tributed rates governed by the distribution of electron-hole radii. Quantum-chemical calculations 
suggest that hot-exciton dissociation at such interfaces results from a high charge-transfer charac- 
ter. 

PACS numbers: 72.20.Jv, 78.47.jd, 78.55.Kz, 78.66.Qn 
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Unravelling primary electronic processes in polymeric semiconductors opens a fundamen- 
tal window to their materials physics, which is of central importance for emerging applica- 
tions in optoelectronics. The steps to generate charge by optical absorption are currently 
the subject of wide interest Here, we focus on charge generation and recombination dy- 
namics in neat regioregular poly(3-hexylthiophene) (P3HT). This semicrystalline polymer 
adopts 7r-stacked lamellar microstructures in the solid state [2], leading to two-dimensional 
electronic dispersion J3]. Crystallinity induced by molecular organization influences pro- 
foundly electronic properties, exemplified by the high yield (77) of apparently direct charge 
photogeneration. Whereas in general rj <C 1% in less organized polymeric semiconduc- 
tors |4| , various groups have reported 77 up to 30% over ultrafast timescales in P3HT films at 
300 K 5|-|9| . This is very surprising because photoemission spectroscopy measurements place 
the energy for free polaron generation ~ 0.7 eV above the optical gap in neat regioregular 



P3HT 



10). Weak interchain electronic coupling in the lamellar architecture leads to a free- 



exciton bandwidth — the pure electronic bandwidth due to dispersion neglecting coupling 



to vibrations — that is well below this energy offset in P3HT stacks 11 



Thus, are 



charges indeed generated intrinsically (directly by photoexcitation, resulting from the crys- 
talline electronic structure) or extrinsically (due to a driving force for dissociation following 
exciton formation)? 

Charge photogeneration in neat P3HT has been studied previously by means of transient 
absorption spectroscopy , which measures directly the dynamics of nascent polarons. 
In this letter, we implement time-resolved photoluminescence (PL) spectroscopy at 10 K to 
probe charge photogeneration and recombination dynamics. Our strategy is to exploit the 
intricate detail of electronic structure, structural relaxation, and correlated disorder afforded 
by the evolution of the spectral bandshape of the PL spectrum [12], over timescales where 
it is known that charge photogeneration is important. By analyzing the bandshape and de- 
cay dynamics of delayed PL due to charge-pair recombination, we probe the recombination 
environment, which also probes the environment in which photogeneration occurs, since pho- 
toexcitations are frozen at 10 K. We find that charge generation occurs over sub-nanosecond 
timescales by dissociation of excitons created at interfaces between lamellar (aggregate) and 
poorly stacked (non-aggregate) domains, driven by energetic disorder. Thus, prompt charge 
photogeneration is an extrinsic process, and film microstructure determines the surface area 
and the energy landscape of interfaces between domains. 
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PL measurements were carried out with a 40-fs, 532-nm (2.33-eV) pulse train derived from 
an optical parametric amplifier (Light Conversion TOPAS), pumped by a Ti:sapphire laser 
system (KMLabs Dragon). Time- resolved PL spectra were measured with an optically- 
triggered streak camera (Axis-Photonique, 6-ps instrument response). Alternatively, an 
intensified CCD camera (Princeton-Instruments PIMAX 1024HB) was used. P3HT films 
(Merck, M w = 47.8kg/mol, M n = 26.2kg/mol, polydispersity = 1.83, 150-nm thick) were 
spun from trichlorobenzene solution (6% wt). 

In P3HT films, the PL spectrum is understood within the framework of a weakly-coupled 
H-aggregate model |l3j . resulting from weak resonance-Coulomb coupling (J) of transition 



moments aligned co-facially in neighboring polymer chains ll). In contrast, the absorption 



spectrum contains contributions from both the aggregate and non-aggregate regions 13]. 
From the ratio of the (0,0) and (0,1) absorbance peaks 14], we estimate a free-exciton 
bandwidth (W = 4 J) of 100 ± 3meV L^. Furthermore, using a Franck-Condon model for 
the absorption spectrum of the H aggregate fl^ . we estimate a volume fraction 52 ± 5% for 
aggregates (electronically coupled chains) versus non-aggregate regions. 

We carried out time-resolved PL measurements at 10 K to explore early-time exciton 
dynamics in neat P3HT, shown in Fig. []Ja). The PL intensity decays substantially over sub- 
nanosecond timescales, and the spectrum red-shifts by > 40 meV with weak evolution of the 
spectral bandshape. We examine in more detail this spectral evolution in Fig.[T^b). Fig. []Jc) 
displays the early-time PL anisotropy, defined as r = — I±)/(I\\ + 2I±), where J||(j_) is the 
instantaneous PL intensity parallel (perpendicular) to the linear polarization of the pump 
pulse. At 10 K (open circles), r decays to 0.23 within the instrument response of several 



picoseconds due to exciton self trapping in the photophysical aggregate [la, Il3|, and does 
not evolve further on a sub-nanosecond window. In contrast, r decays to zero after 500 ps at 
room temperature. PL anisotropy decay is a signature of exciton diffusion [19], so at 10 K, 
excitons are immobile over the timescale of the dynamic red shift in Fig. [T^a). Another 
possibility could be relaxation of excess vibrational energy following ultrafast excitation 
of the sample. Parkinson et al. have reported that torsional relaxation of the backbone 
leads to more cofacial lamellar structures (more extended correlation of site energies in the 



aggregate [12]), leading to a dynamic loss of (0,0) intensity over ~ 13 ps at 300 K [l8[. We do 
not observe significant relaxation of the (0,0) relative intensity over the timescale of Fig[]Ja). 
Furthermore, upon increasing the temperature even slightly, the steady-state PL spectrum 




FIG. 1. (Color online) (a) Time-resolved PL spectrum, exciting at 2.33 eV with fluence 19 ^Jcm -2 , 
measured at 10 K. (b) Normalized slices of (a) at times indicated in the legend, (c) Time-resolved 
PL anisotropy r at 10 K (open circles) and at room temperature (dots). The horizontal line 
indicates the average anisotropy, r = 0.22 ± 0.03, between 20 and 600 ps. 

broadens significantly at 10 K [13J, so we rule out thermal relaxation. We suggest that 
the dynamic spectral red-shift of the PL spectrum results from an evolving electric field 
distribution [20j from photogenerated charges in the aggregate, which is consistent with 
charge photogeneration yields reported in the literature and our own transient absorption 



measurements 14| • If so, charges are generated over all timescales spanning < 1 ns, by a 
mechanism not involving exciton diffusion. 

We next consider delayed PL dynamics due to charge recombination. Fig. [2]^a) displays 
time-gated PL spectra at 10 K. After a few nanoseconds, the aggregate spectrum is super- 
imposed with a broader component (see ref [ljj]) with similar dynamics to those of the H 
aggregate since this composite bandshape persists over microseconds. While we know that 




FIG. 2. (Color online) (a) Time-resolved PL spectra measured after 2.33-eV excitation at 10 K, and 
at gate times denoted in the legend, (b) Spectrally integrated PL intensity versus gate delay. The 
time-resolved PL intensity from Fig. [pjis included over sub-nanosecond timescales. The dashed lines 
represent a two-exponential fit with (r) = 235±23 ps and an asympototic power law [l/(H-(t/r / ) Q )] 
with a = 1.54 ± 0.08 and t' <C (t). A fit to the model (eq. [1]), represented schematically in 
(c), is shown as a continuous line (for the case of r\ = 0.03, T c t = 0.35 ns -1 , F x = 1.23 ns _1 , 
v = 4.37 x 10 13 s , ijl = 0.54). (d) Delayed PL excitationspectrum at 10 K, measured by phase- 



sensitive methods as described in auxiliary information [14|, superimposed on the components of 
the H-aggregates and non-aggregates of the absorption spectrum (see ref. Il5l). Also shown is the 
external quantum efficiency spectrum measured at 300 K in a ITO/PEDOT:PSS/P3HT/LiF/Al 
diode fabricated and encapsulated in a N2-filled glove box. 



the aggregate spectrum arises from lamellar stacks, we conjecture that the featureless spec- 
trum arises from non-aggregate regions. It is reminiscent of red-emitting species found in 
amorphous polymer films, often referred to as excimers 21] . From the delayed PL bandshape 
in Fig. |5£a), we conclude that recombination occurs at the interface between H-aggregate 
and non-aggregate regions, and can populate either species. 

The total PL intensity decays exponentially over picosecond timescales, but asymptoti- 
cally as a power law (I(t) oc t~ L54 ) for times much longer than the exciton lifetime (Fig.|2fb)). 
This behavior can arise 
ular annihilation 

pump fluence [14j, which permits us to assign the power-law decay to exciton regeneration 



23 



rom unimolecular charge recombination [22j or by triplet bimolec 



241 ] . We find that the delayed PL intensity varies linearly with the 



by charge recombination with a distribution of rates, as the other two possibilities would 
lead to a non-linear fluence dependence. The power-law decay is independent of tempera- 
ture 14J, suggesting recombination by tunnelling. We note that we cannot reproduce the 
measured time-resolved PL intensity as a simple superposition of a multiexponential and 
an asymptotic power-law decay, as the time window spanning 1-10 ns features a PL decay 
over nearly three orders of magnitude that deviates significantly from either decay function 
(Fig. Hl^b))- This indicates that the two PL decay phenomena are not independent, but that 
one decay regime evolves into the other, with competing kinetics on nanosecond timescales. 
By integrating the PL intensity over timescales where the decay is non-exponential, we find 
that > 12% of the time- integrated intensity is accounted for by slow recombination. This 
reflects a significant density of charge pairs. 



Based on the results of Figs. [T] and [2] and building upon previous reports 



we con- 



struct the following photophysical picture, depicted schematically in Fig. Upon pho- 

toexcitation, charge pairs are generated with efficiency i], and the rest of the population 
relaxes to self-trapped excitons x in <lps [l^]. These decay to the ground state with rate 
constant T x , or charge-separate with rate constant T ct . Charge pairs may recombine to 
regenerate x with a temporal distribution R(t). Thus, 



d.r 
dt 



-(T x + T ct )x + V R(t)+ [ T ct x(t')R(t-t')dt'. (1) 

Jo 



Assume a distribution of charge-pair distances /(r) = ee~ er . The tunnelling-rate distance de- 
pendence is k{r) = z/e _/3r and the tunnelling-time distribution is R(t) = J °° f(r)k(r)e~ k ^ r ' t dr. 
With x(0) — 1 — 7], the Laplace transform of eq. [T] is 

m = — T , ( 2 ) 



r. + r 



ct 



R(s) 



where R(s) = \i J °° e _<1+At ) x /(s/V + e~ x )dx, X — P r i an d ^ = The model predicts 

biphasic dynamics with PL intensity, I(t) oc T x x(t), evolving as I(t) oc exp(— T x t) at short 
times and as I(t) oc t~^ 1+ ^ at long times. We evaluate eq. [2] as described elsewhere to 
recover I(t) [2q . 

The most robust parameter in this model is the ratio of the characteristic electron-hole 
separation (e) and the distance dependence of charge tunnelling /i = e/(3 = 0.54 ±0.08, 
as it defines uniquely the slope of the power-law decay at long times. Hence, e ~ 0.54/3; a 



significant fraction of charge pairs do not recombine in the microsecond timescale, as the tail 
of f(r) extends beyond the characteristic lengthscale of k(r). The charge-pair population 



is expected to survive on much longer timescales, w 



rich is consistent with reports of a high 



26 



27). 



residual charge density in P3HT at steady state 

We next consider 77. We cannot extract it uniquely from the model, as it is coupled 
to T ct . Both parameters affect the amplitude of the power- law decay without altering the 
decay rate. We probe the range of rj that is consistent with the literature, and explore the 
consequences on T ct . With = 1.23 ns~ x (fixed by the slow part of the bi-exponential decay 
in Fig. Ufa)) and with v = 4.37 x 10 13 s _1 (corresponding to the frequency of the aromatic 
C — C stretch measured by Raman spectroscopy [28]), we find that if we set rj = 3%, we 
require T ct = 0.35 ± 0.05 ns _1 , but if rj = 30% we can reproduce the data with T ct = 0. 
With 77 = 40%, we can no longer obtain a satisfactory fit of the amplitude of the power- 
law component, which is over-estimated. We have interpreted the dynamic red shift of the 
PL spectrum (Fig. [TJ as a consequence of an evolving electric field on the subnanosecond 
timescale, which is comparable to the exciton lifetime. In order to rationalize our time- 
resolved spectroscopic data, we therefore consider that T ct ?S T x , which from the modelling, 
would imply that r\ < 10% in the solid-state micro structure of our films. 

If the two slow-decaying emissive species in Fig. E(a) are due to regeneration of aggregate 
and non-aggregate excitations by charge tunnelling, we propose that charges are initially 
generated at the interface between the two domains. To explore this further, we have carried 
out quantum chemical calculations j^J on a stack of ten oligothiophene chains (Fig. 13(a)) 
to represent the crystalline moieties in the higher molecular-weight (i.e. longer chain) P3HT 
used in our experimental study. We note higher conformational disorder at the top/bottom 
of the stack, leading to a higher excitation energy in those regions compared to the center 
(Fig. E(b)). Fig. [3(c) displays the calculated absorption spectrum in a given configuration. 
For all stack configurations studied, we always observe a shoulder at ~ 2.7 eV on the blue side 
of the main absorption band at ~ 2.4 eV. (Note that these energies are overestimated, but 
the ~ 0.3 eV difference is meaningful.) The excited states that are generated in this spectral 
region are quasi-degenerate with the lowest-lying charge-transfer (CT) states, shown by the 
superimposed plot of CT character as a function of excitation energy. We can identify three 
regimes of the spatial distribution of transition densities, shown in Fig. [3(d). The excitation 
with the lowest transition energy (I) is always confined to the center of the stack over two 



FIG. 3. (Color online) Results of quantum-chemical calculations [29]. (a) Snapshot of the archi- 
tecture of the aggregate (15-mers in a stack of 10 chains), (b) Average excitation energy (error 
bar is the variance) for the (10) chains in the stack, averaged over a large number of excited-state 
calculations on isolated conformers extracted from the stack, (c) Absorption spectrum for one 
configuration of the full stack superimposed to the corresponding charge-transfer character, (d) 
Transition density diagrams (ground-state - excited-state overlap) for three excitations of increas- 
ing energy. 

to three sites as a result of disorder and is weakly emissive (H aggregate); the intermediate 
excitation (II) carries most of the oscillator strength and is delocalized in the center of the 
stack; the higher-lying excitation (III) is dominated by conformationally disordered chains at 
the ends of the stack, but it communicates with chains in the center so that charge separation 
of these higher-energy states is possible. In our measurements, we excite ~ 0.3 eV above 
the (0,0) absorption, placing us at the limit of CT excitation predicted by Fig. [3]^c). We 
have measured the delayed PL excitation spectrum, shown in Fig[2]^d). This reveals that 
excitation on the edge of the component of the absorption spectrum due to less ordered 
or amorphous species enhances the delayed PL yield, consistent with the significant CT 
character of excitons at these photon energies. Photocurrent measurements reveal that a 
further 0.4-eV energy is required to produce charge carriers (Fig. EJ^d)), underlining that our 
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delayed PL measurements probe recombination of tightly-bound geminate polaron pairs (not 
carriers) produced by excitation with energies ^ 3eV. Deibel et al. reported that 0.42-eV 
excess energy above the optical gap is required to produce polaron pairs, and then a further 



0.3 eV to generate photocarriers 



101 ]. consistent with our findings. 



We have presented compelling evidence that excitons with sufficient energy dissociate 
extrinsically near interfaces between molecularly organized and less organized domains, un- 
derlining the key role of the disordered energy landscape in semicrystalline microstructures. 
This points to the importance of hot-exciton dissociation in neat P3HT, which competes 
with dissociation at donor-acceptor heterojunctions in photovoltaic diodes [l|. 
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